Stability of the two-dimensional Bose gases in the resonant regime 
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We consider a two-dimensional Bose gas formed in a planar atomic trap in conditions where the 
two-dimensional scattering length exceeds all other microscopic length scales in the system and, 
accordingly, the gas parameter assumes relatively high values. We show that, unlike in the three- 
dimensional case, for sufficiently low areal densities the two-dimensional gas remains stable against 
collapse even in the resonant regime. Furthermore, we evaluate the three-body recombination rate 
that sets the upper limit for the life-time of two-dimensional resonant atomic condensate. 
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I. INTRODUCTION 

Since the achievement of Bose Einstein Condensation 
(BEC) of neutral atoms in 1995, the field of cold atoms 
has been substantially developed, both in the parameter 
range (number of atoms, temperature, etc), and in the va- 
riety of systems studied. In particular, low-dimensional 
atomic systems became a subject of significant interest 
as a consequence of the existence of nontrivial strong 
phase fluctuations and also due to the peculiar scatter- 
ing properties Two-dimensional (2D) Bose gases play 
a special role in this sub-field thanks to a rich phase 
diagram Q. Experimentaly, a two-dimensional quasi- 
condensate has been first achieved with spin-polarized 
hydrogen atoms adsorbed on a liquid 4 He film [3[. Sev- 
eral groups have now reached the 2D regime with ultra 
cold atoms trapped using evanescent light waves [ij or 
optical lattices |5j. In the latter trap, thermaly activated 
vortices which plays a central role in the Berezinsky- 
Kosterlitz-Thouless phase transition have been observed 
experimentaly Q. Ultra-cold atoms in two dimensions 
are also interesting due to the fact that the atom-atom 
effective coupling constant can be easily tuned 0, Q by 
varying the axial confinement of the trap and/or by us- 
ing a Feshbach resonance Consequently, ultra-cold 
atoms experiments open interesting perspectives for a 
detailed study of atomic BEC in strictly (i.e. with no 
excited transverse modes present) two-dimensional con- 
figurations in a broad range of experimental parameters. 

Using the 2D equation of state, we show that in the 
resonant regime where the two-dimensional scattering 
length exceeds the size of the transverse confinement 
(and becomes the only relevant length scale associated 
with the interatomic interactions) the gas remains sta- 
ble (at least for sufficiently small areal densities), even 
though the corresponding three-dimensional scattering 
length becomes negative in this regime. Moreover, since 
in the resonant regime the 2D gas parameter becomes rel- 
atively high, the formation of dimers in three-body colli- 
sions becomes relevant |l2j. Therefore, we study further 
the stability of the resonant two-dimensional condensate 
vis a vis this process. 



The paper is composed of two parts. In the first part, 
we review common low energy collisional features for two 
ultra-cold atoms in a planar harmonic guide and give the 
definition of the 2D scattering length (a^rj)- The effective 
interaction between two confined atoms is described us- 
ing the zero-range approximation within the A-potential 
formalism. The resonant regime for two-body scattering 
is defined by the condition a2_D 3> a z . A bound state of 
vanishing energy and large spatial extension (of the order 
of a 2 c) emerges under these conditions: this is nothing 
else but a dimer which is populated by three body pro- 
cesses in a dilute resonant gas. 

In the second part, we consider the 2D atomic Bose 
gas of areal density n. For low areal density systems 
where the interparticle distance is much larger than the 
transverse confinement na 2 , <C 1, the interaction can be 
described in terms of the 2D scattering length only. For 
small gas parameters na^ D <C 1, the two-dimensional 
equation of state shows no sign of a collapse even though 
the corresponding three-dimensional scattering length 
azD is negative in the two-dimensional resonant regime. 
Finally, we evaluate the rate of three-body atomic recom- 
binaison into the shalow dimer and give an estimate for 
the life time of a two-dimensional atomic condensate at 
the onset of the resonant regime. 



II. TWO-BODY SCATTERING IN A 
TWO-DIMENSIONAL PLANAR WAVE-GUIDE 

A. Planar wave-guide 

We consider two atoms in an harmonic planar waveg- 
uide characterized by the angular atomic frequency lo z 
in the transverse direction z. The characteristic length 
is given by the size of the ground transverse mode a z = 

— . The two atoms have the same mass m and are 
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located at positions Ri and R2, respectively. We further 
decompose the position vectors into a sum of longitudinal 
(parallel to xy-plane) and transverse (along z axis) com- 
ponents: Ri = ri -\-z\e z and R2 = Y2 + Z2e z , where e z is 
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a unit vector in the transverse direction. This system can 
be considered as 2D when the energy of the relative mo- 
tion E is close to the threshold between the monomode 
and non-propagating regimes: e = E — hw z /2 <C hio z . 
Then, for |ri — r 2 | ^> a z , the two-body wave function 
factorizes as: 

* 2 (Ri,R 2 ) =00(^1)^0(^2)^2^1^2), (1) 



B. Two-dimensional scattering amplitude and 
two-dimensional scattering length 

In a scattering process at energy e = h 2 k 2 /m (with 
ka z <C 1), and for large relative distance (r 3> a z ), the 
2D wave function is 



1P2 ( r ) = exp(ik.r) 



2 ln(— ikqa 2 o) 



ijW(fcr). (2) 



In Eq.@ r = ri — r 2 are the relative coordinates, the 



constant a 2 u is the "2D scattering length" and ( 



A 



where 7 is the Euler constant. It can be shown [8|,0 that 
a 2D is a function of the characteristic length of the trap 
a z and also of the three dimensional scattering length 
a 3D with [l(J: 



a>2D = 2.092 a z exp 



\ ±- ) • (3) 

2 a 3D 



C. Zero-range approximation in two dimensions 

Let us introduce the notion of a zero-range two- 
dimensional interaction potential for the two trapped 
atoms. For this purpose, we use the fact that the radius 
of the interatomic forces R is very small as compared 
to all the lengths considered in this paper, so that the 
3D wave function is a solution of the free Schrodinger 
equation (r ^> R): 



ti 2 1 
— (Ai + A 2 )* 2 + -mcj 2 z {z 2 + 4)* 2 



E^n 



(4) 



Consequently, for r 3> a z and e <C ftiu z , the 2D wave func- 
tion ip2 is also solution of the free Schrodinger equation. 
Hence, "02 has a logarithmic behavior for small interpar- 
ticle distances (in the scattering process this corresponds 
to fcr C 1 but r ^> a z ). The idea of the zero-range 
approach is to consider the extension of ■02 defined as 
a solution of the free Schrodinger equation Vr > and 
to impose the correct behavior of the wave function for 
r»Oj from a contact condition taken formally at r = 0. 
Using the expression of the scattering states in Eq.([2|) 
extended in the region kr <C 1, one can verify that all 
the two-body states which are linear combinations of $Z$i 
can be deduced from the contact condition: 



■02(i"i,r 2 ) = A 2 \n[ 

V CL2D 



■0(r), 



(5) 



where A 2 depends on the state considered and also can 
be a function of the center of mass coordinates. The zero 
range two-dimensional potential introduced in Ref.(l3l] 
allows to incorporate the contact condition ([5]) in the 
Schrodinger equation defined Vr > (at r — 0, the iden- 
tity Aln(r) = 27T(^ 2 )(r) is used). It has the following 
expression: 



2ttH 2 S (2 \r) 



log(Aa 2r j) 



1 — log(Ar)r 



0_ 

dr 



, (6) 



where A is an arbitrary constant. This potential acts on 
wave functions of the form ([5|) as: 



(ri,r 2 |V$|V*> 



(7) 



and gives a result independent of A. Let us emphasize 
that this property is not general: the action of on 
a two-body wave function which has not the singularity 
given by Eq.© leads to a A-dependent result, showing 
also that such wave function is not in the proper Hilbert 
space. 



D. Resonant regime: shallow bound state 

In the regime of large 2D scattering length a 2 D 3* a z 
which corresponds to the situation of negative 3D scat- 
tering length dsn < 0, the scattering cross section is max- 
imum for k — kb = (qci 2 d) _1 . In this resonant regime, 
there exists a low energy two-body bound state (refered 
hereafter as the dimer state) with a spatial extension of 
the order of cl 2 d- Its wave function (0b) and binding- 
energy (es) can be deduced from the contact condition 
([5]) or also from the pseudo-potential ([6]). One obtains: 



f>B(n.2) 



KB 



K (K B ri 2 ) , and e B 



h 2 n 2 B 



(8) 



Now, the key point for the following part is that the 
state in Eq.© has an extension of the order of a 2 u- as 
a consequence it corresponds to the asymptotic part of a 
physical one in the regime: 



a 2 D > a z 



(9) 



As a consequence of the exponential behavior in Eq.Q, 
this resonant regime is achieved for negative three dimen- 
sional scattering length and: 



< -a 3D < a z . 



(10) 



III. IS THE DILUTE RESONANT 
TWO-DIMENSIONAL BOSE GAS STABLE? 

A. Equation of state: stability against collapse 

We consider an atomic Bose gas characterized by an 
areal density n and confined in the planar wave guide 
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with na 2 <C 1, so that only the ground transverse mode 
is populated. In the dilute limit defined by a small 2D 
gas parameter: 

na^<l, (11) 
the 2D equation of state at zero temperature is [3, HE 

mfj, ( 4h 2 \ 

H= 4^ ( ve (27+D mAla 2 D J ■ ( 12 > 

For sufficiently small areal densities such that Eq. fTTj) 
is satisfied, Eq. (TT2")) shows that mpa 2 D /h 2 <C 1 and 
dn/dfi > 0, Va2£). Consequently for a vanishing areal 
density, the system is thermodynamically stable which 
means an absence of collapse. This is in particular 
true in the regime of negative three dimensional scat- 
tering length a^D < 0, where an instability occurs for 
the homogeneous gas in absence of trapping. Due to 
the exponential behavior of a 2 £> with respect to a 3 £> in 
Eq.Q, the condition (fTTj) is too stringent to reach for 
< — <Z3£) <C a z and the more interesting situation cor- 
responds to the onset of the resonant regime: 

< -a 3D < a z . (13) 

However one can notice that even at the onset, the condi- 
tion defined by Eq. (fTTj) is severe and is difficult to satisfy 
for clouds involving thousand of atoms (as a consequence 
of the radial trapping) [llj • 

However, even if the 2D Bose homogeneous gas is sta- 
ble against collapse in the peculiar regime defined by 
Eqs. pill3p . its stability is limited by three body recom- 
binaison processes into the shallow dimer Eq.©. In the 
next section we adress this question and give an estima- 
tion for the life-time in this regime. Hence, we give a 
limit on the possible occurence of an atomic condensate 
for increasing values of the 2D gas parameter obtained 
by a variation of the 2D scattering length. 

B. Three-body recombination rate: stability 
against formation of dimers 

The problem that we are considering is rcmincsccnt 
of the three dimensionnal resonant Bose gas where three 
body recombinations populate a dimer state of vanish- 
ing energy and diverging width [l^]. More precisely, the 
transverse length a z plays the role of an effective range 
and Eq.© is analogous to the resonant regime in 3D 
where a^D ^ R. However the two situations are not com- 
pletely similar. Indeed, in the 2D resonant regime Eq.© 
a single parameter a 2 £> is sufficient for a characterization 
of the low energy properties in few- and many-body sys- 
tems [19f. while in 3D at least one more parameter is 
needed [201 - Moreover, in 2D the low energy two-body 
cross section strongly depends on the energy so that one 
expects a strong variation of the recombination rate with 



respect to the gas parameter which fixes the typical col- 
lisional energy. A straighforward dimensional analysis 
gives the following law for the recombination constant: 

a rec = F{na\ D )—a 2 2D , (14) 
m 

and the purpose of the subsequent three-body analysis is 
to evaluate the function F. 

We consider three atoms (1,2,3) of equal mass m 
which after a collision lead to the formation of a dimer 
[2lj |. There are three outgoing channel for this process 
obtained by cyclic permutation from: 

1 + 2 + 3^1 + (23) (C), (15) 

where (23) denotes the dimer and (C) labels this par- 
ticular channel. The problem can be solved using the 
Faddeev equations [22[ in spatial coordinates. This ap- 
proach is very efficient for the determination in the con- 
figuration space of the two universal three-body bound 
states in two dimensions |23j (these states where first ob- 
tained in the momentum representation in Ref.[24]). The 
positions of the atoms (1, 2, 3) are defined on the plane by 
the coordinates (ri,r 2 ,r 3 ) respectively. The wave func- 
tion factorizes in a z-part frozen in the ground state of 
the transverse harmonic oscillator (omitted in the follow- 
ing) and a 2D part 5' 3 (r 1 , r 2 , r 3 ). The wave function $3 
is an eigenstate of the Hamiltonian: 

^3 = - E + V 1 A 2 + V li + V 2% (16) 

To be consistent with the 2D regime, the eigen-energy 
of this state must be below the transverse level spacing, 
e <C huj z . The potential part in Eq. (|16[) imposes the con- 
tact condition for cyclic permutations of the three 
particles. In the center of mass frame, the wave function, 
which is symmetric as a consequence of the symmetriza- 
tion postulate is decomposed as: 

¥ 3 (1, 2, 3) = r/>(ui, u 2 ) + V(ui, u 2 ) + u 2 ') , (17) 

with the Jacobi variables defined by: 

"l = \/f( r i ~ ~ 2 ~ ) ' u 2 = ^|( r 2 ~ r 3)- (18) 

The other variables (u' 1; u 2 ) and (u", u 2 ) can be also in- 
troduced via cyclic permutations of the three atoms. We 
use then the hyperspherical coordinates {p, a, a', a"} de- 
fined by: 

Ui = pcos(a) ; u 2 = psin(a) , (19) 

and hyperangles a 1 and a" are obtained by cyclic permu- 
tations. As we are interested in a low energy process, the 
3-body wave function is a s-state which does not depend 
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on the planar angles. Consequently, the kinetic operator 
in Eq. p6p is simply 



T = ^\-^d 2 ( P ^ 2 .) + ^^ 



2m | p3/2 P 



P* 



with, A 2 = -d 2 



(20) 



(21) 



tan(2a) 

In this coordinate system, the contact condition reads 
* 3 = A 3 (p)]n(^-)+G(a), (22) 

y 0.2D J 

and the contact conditions for a' and a" have an analo- 
gous form. We are going to show now, that the knowledge 
of the function A 3 (p) is sufficient for the determination 
of the three body recombination rate. For this purpose, 
we use the expression of the on-shell i-matrix element 
from the initial channel of three atoms in the continuum 
of positive total energy e, to the outgoing channel (C) 



*e = <*(C)|V(c)|*3> 



(23) 



In Eq . (|2"3"|) . $3 is the exact 3-body eigen- function of the 
Hamiltonian (|16|) corresponding to the initial state, and 
* s ^ ne non-interacting 3-body wave function for the 
outgoing channel (C): 



Ui,u 2 |#$) = — exp ( ik.un/- ) <M|r 2 -r 3 |), 



V3 



where the wave number k is obtained from the energy 

conservation in the 3-body process: k — — ■=■. The in- 

v3 

teracting potential (V(c)) which appears in Eq. (j2"3")) is the 
scattering potential for channel (C) [25| : 



V { c) = V& + V& . 



(25) 



Summing up the contributions of the three possible out- 
going channels, one obtains the total t-matrix element: 

*i/ = 6<*$|^|*3>. (26) 



From Eqs. (|22l26[) . the total t-matrix element can be de- 
duced from: 



if 



4ir 2 V3h 2 



dp pJa{-^)<t>B{p 



)MP)- (27) 



This relation can be understood qualitatively if one re- 
marks from Eq. ([22|) . that A 3 (p) is linked to the probabil- 
ity of finding one atom at distance p from a colliding pair 
with an interparticle spacing smaller than cl2d- Then, as 
the integration of this function in (|27|) is over lengths 
of the order of a 2 £>, the value of tif depends crucially 
on the probability of finding three atoms at relative dis- 
tance 02D, which coincides also with the extension of the 



dimer. Note that the result (f27| does not depend on A as 
expected for a theory which provides the exact contact 
conditions on the wave function. 

The determination of the function A3 can be in princi- 
ple obtained from the decomposition of the wave function 
over a hyperspherical harmonics basis [23| : 



*3 = E 



Mp) 



p 



,3/2 



$i(jO, a, a', a") 



(28) 



where a, a', a") = <pi(p,a) + 4>i(p,a') + <fii(p,a") 

and for each value of the hyper-radius p, tfii (p, a) is an 
eigen-function of the operator A 2 with the eigen-value 
\i(p) (this step is very reminiscent of the 2-body prob- 
lem). Taking into account the fact that the functions 

7T 

(f>i (p, a) are finite for a — — , we obtain: 

4>i(p,a) =AfiP Vi (-cos(2a)), (29) 
where V Vi is the usual Legendre function and: 



-i + yr+Mp) 



(30) 



For a given value of p, the normalization constant Mi 
is defined by integration over the angles a and (6*1,6*2) 
(associated to (ui,u 2 )): 

\ da d6 1 I dd 2 |$i| 2 sin(a)cos(a) = 1. (31) 
Jo Jo Jo 



In the limit a — > 0, the functions </>j in Eq. (l29|) exhibit a 
logarithmic singularity: 

^(p,a) = Msin(7ri/i) ( ln(a 2 ) + 2 7 + 2^ i + l) 



+ 0(a), (32) 



tan(7r^:) 



where ^(v) is the logarithmic derivative of the gamma 
function. Hence, for each value of p the contact condition 
in Eo.(|2"2"|) is satified if the eigenvalue \{p) is a solution 
of the equation: 



j + ip(vi + 1) + 



2tan(-7Wi) sin(7rz/i) 



^(l/2) =l JpV2\ 33) 



0-2D 



The left hand side of Eq. (f3"3")l is a function of Xi that we 
denote by L(Xi). This function has several branches asso- 
ciated to the harmonics in Eq. (|28|) and the two first ones 
are represented in Fig.fl]). Solution of Eq. (j3U|) in each 
branch defines a hyperspherical harmonic $i. As the ini- 
tial non interacting state corresponds to three atoms in 
the continuum with a vanishing energy, the most impor- 
tant weight in the decomposition of Eq. (|2"5|) is the har- 
monic corresponding to the lowest positive eigenvalue. 
We then limit our study to a single harmonic analysis 
with the wave function: 



n 3/2 



4>(p, a, a' , a") . 



(34) 
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FIG. 1: The two lowest branches of the function L - left hand 
side of Eg.(|33jl. 



The approximation used here is linked to the fact that 
we can identify a slow variable (p) and a fast variable 
(a) which are almost separable. The function A^{p) can 
then be deduced from: 



2/(p) 



3/2 



\(ttv)M{p) . 



(35) 



For each value of p > 0, we solve Eg . (|3"3")) in the do- 
main A £ [0,24]. This way we construct the function 
X(p) which leads to an effective kinetic barrier in the one 
dimensional equation for f(p) : 



f = K 2 f . 



(36) 



In Eq. ([36|) . K is the momentum defined by the energy 
of the process e = h „ K (Ka,2D C 1). In the limit of 



very large hyper-radius, the amplitude of / is fixed by 
considering the asymptotic behavior of the free incoming 
wavefunction, so that 



/ 67T 

p— >oo y i\ ,J 



(37) 



We apply now this formalism by considering three atoms 
initially in the 2D atomic Bose condensate of areal den- 
sity n at the onset of the resonant regime (fT3)) . The 
typical energy in the collision process is given by e ~ 3p. 
For simplicity the chemical pot ential can be evaluated by 
use of the Schick's formula [26j which is an approximate 
form of the equation of state: 



A-n:h 2 n 



p 



m\\n(nal D )\ 



(38) 



The gas parameter is small na\ D <C 1 which implies also 
that the collisional energy is such that e -C \es\ <C huj z . 
An analysis similar to Ref.[l8| gives the number of pro- 
cesses per unit time and surface and permits an iden- 
tification of the function F(x) (where x = na\ D is the 
2D gas parameter) entering in the recombination con- 
stant a r 



CU with F{x) = m 2 tj f l{m 2 a 2 D) 2 



and the 
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2d 

FIG. 2: Function l/(x 2 F(x)) entering in Eq.((39j 

life-time of the atomic Bose gas is : 
t{x) 



ma 



2D 



3hx 2 F(x) 



(39) 



In Fig.@, we have plot the function l/(x 2 F(x)). 

As an example, we consider 7 Li atoms confined in 
a wave-guide of atomic frequency lo z = 2it x 250 kHz 
(a z — 76 nm and is very much larger than the range 
of the true interatomic forces ~ 2.8 nm) with a 2D gas 



parameter x 



L 2D 



= 5x 10 2 cm 2 and a 2 r> = 10a z . 



We obtain for such configuration a very short life time 
r ~ 2 x 10~ 2 s, showing that the 2D atomic condensate 
decays rapidly by creating dimers in three-body 2D re- 
combination processes and not by experiencing a collapse 
phenomenon (a^D < 0). For a smaller but also "large" 
2D gas parameter as x = 10~ 2 with a2D = 10a z , then 
r ~ 3 s showing that the atomic condenstae may be ob- 
servable for relatively high values of the gas parameter. 
However, as shown on Fig. J2]), l/(x 2 F) exhibits a rapid 
decrease with respect to the gas parameter x, so that for 
x on the order of 0.1, the atomic condensate decays very 
rapidly. 



IV. CONCLUSION 

In this paper we analyse two potential sources of insta- 
bility of two-dimensional atomic Bose condensates in the 
universal resonant regime of the large two-dimensional 
scattering length (corresponding to a negative three- 
dimensional scatterring length of the order or less than 
the size of the transverse confinement): a macroscopic 
collapse and a loss of atoms through formation of shal- 
low dimers in three-body collisions. We show that (a) the 
gas is stable against collapse and that (b) in the above 
regime the upper limit on the condensate life-time set 
by the three-body recombination can be relatively long. 
Moreover, for a trap sufficiently deep in the longitudi- 
nal direction, dimers can be kept in the trap opening the 
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possibility of a series of recombination process with suc- 
cessive formations of three- and four-body bound states 
toward "2D Bosons droplets" states such as predicted in 
Ref.[I3. 
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